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Abstract 

We consider the multivariate interlace polynomial introduced by Courcelle 
(2008), which generalizes several interlace polynomials defined by Arratia, Bol- 
lobas, and Sorkin (2004) and by Aigner and van der Hoist (2004). We present 
an algorithm to evaluate the multivariate interlace polynomial of a graph with 
n vertices given a tree decomposition of the graph of width k. The best previ- 
ously known result (Courcelle 2008) employs a general logical framework and 
leads to an algorithm with running time f{k) ■ n, where f{k) is doubly expo- 
nential in k. Analyzing the GF(2)-rank of adjacency matrices in the context of 
tree decompositions, we give a faster and more direct algorithm. Our algorithm 
uses 2^^ . arithmetic operations and can be efficiently implemented in 

parallel. 

1 Introduction 

Inspired by some counting problem arising from DNA sequencing [ABCSOO] , Arratia, 
BoUobas, and Sorkin defined a graph polynomial which they called interlace poly- 
nomial |ABS04a] . It turned out that the interlace polynomial is related |ABS04a| 
Theorem 24] to the Martin polynomial, which counts the number of edge partitions 
of a graph into circuits. This polynomial has been defined in Martin's thesis from 
1977 |Mar77j and generalized by Las Vergnas |LV83] . Further work on the Martin 
polynomial has been pursued p7M [LV881 1.Tae88 l. iEM98i 'EM99', 'Boin2]. including 
a generalization to isotropic systems |Bou87t [Bou8 8. Bou91, BBD97J. In particular, 
the Tutte polynomial of a planar graph and the Martin polynomial of its medial 
graph are related. This implies a connection between the Tutte polynomial and the 
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interlace polynomial (see Ellis- Monaghan and Sarmiento |EMS07] for an explana- 
tion) . 

One way to define the interlace polynomial is by a recursion that uses a graph 
operation. Arratia et al. used a pivot operation for edges |ABS04aj . This operation 
is a composition of local complementations to neighbor vertices (see Aigner and 
van der Hoist |AvdH04] . where the operations are called switch operations). The 
orbits of graphs under local complementation are related to error- correcting codes 
and quantum states, and so is the interlace polynomial as well |DP08] . 

The interlace polynomial can also be defined by a closed expression using the 
G'-F(2)-rank of adjacency matrices |AvdH04"| IBouOSt IEMS06] . This linear algebra 
approach has been used in several generalizations of the interlace polynomial. In 
this paper, we consider the multivariate interlace polynomial C{G) defined by Cour- 
celle |Cou08] (see Definition 12.11 below) as it subsumes the two-variable interlace 
polynomial of Arratia, Bollobas, and Sorkin |ABS04b] and the weighted versions of 
Traldi |TralO] . as well as the interlace polynomials defined by Aigner and van der 
Hoist |AvdH04] . Furthermore, the interlace polynomials Q{x,y) and Q^^, which 
have emerged from a spectral view on the interlace poljTiomials jRPOGj . are also 
special cases of Courcelle's multivariate interlace polynomial. 

1.1 Results and Related Work 

Our aim is to present an algorithm that, given a graph G = {V, E) and an evalua- 
tion point, i.e. a tuple of numbers {{xa)aeVi {ya)a<^v,u,v), evaluates the multivariate 
interlace polynomial C{G) at {{xa)aev, {yajaeVjUjV). Whereas this is a #P-hard 
problem in general |BH08] . it is fixed parameter tractable with cliquewidth as pa- 
rameter |Cou08[ Theorem 23, Corollary 33]. This is a consequence of the fact that 
the interlace polynomial is monadic second order logic definable {MSi definable as 
defined by Courcelle, Makowsky, and Rotics |CMR01] : see also Courcelle |Cou08t 
Section 5])ll] Such graph polynomials can be evaluated in time f{k) ■ n, where n 
is the number of vertices of the graph and k is the cliquewidth. The function f{k) 

^Note the following crucial difference with respect to monadic second order logic definability: 
MSi definable evaluation problems are fixed parameter tractable with cliquewidth as parameter 
|CMROll Theorem 31]. MSi is a logic that allows one-sorted structures, the universe of which 
consists of the vertices of the graph. Set variables range over vertex subsets. On the contrary, MS2 
is a logic that allows two-sorted structures, the universe of which consists of the vertices and edges 
of the graph. Set variables range over vertex subsets or edge subsets, which, for instance, enables 
the definition of the Tutte polynomial in MS2- MS2 definable evaluation problems are known to be 
fixed parameter tractable with tree-width as parameter ^CMROlt Theorem 32]. We can not expect 
that this generalizes to cliquewidth, see Fomin, Golovach, Lokshtanov, and Saurabh [FGLSIO] . 
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can be very large and is not explicitly stated in most cases. In general, it grows as 
fast as a tower of exponentials the height of which is proportional to the number of 
quantifier alternations in the formula describing the graph polynomial |Cou08l Page 
34]. In the case of the interlace polynomial, this formula involves two quantifier al- 
ternations |Cou08t Lemma 24], |CiO07] . If a graph has tree width k, its cliquewidth 
is bounded by 2^~^^ [COOP] . Thus, the machinery of monadic second order logic 
implies the existence of an algorithm that evaluates the interlace polynomial of an 
n-vertex graph in time f{k) ■ n, where k is the tree width of the graph and f{k) is 
at least doubly exponential in k. (In particular, the interlace polynomial of graphs 
of treewidth 1, that is, of trees, can be evaluated in polynomial time, which also has 
been observed by Traldi |TralO] .) 

The monadic second order logic approach is very general and can be applied not 
only to the interlace polynomial but to a much wider class of graph polynomials 
|CMROlj . However, it does not consider characteristic properties of the actual graph 
polynomial. In this paper, we restrict ourselves to the interlace polynomial so as to 
exploit its specific properties and to gain a more efficient algorithm (Algorithm [2]) . 
Our algorithm performs 2^^^ +o{k)^ arithmetic operations to evaluate Courcelle's mul- 
tivariate interlace polynomial (and thus any other version of the interlace polynomial 
mentioned above) on an n-vertex graph given a tree decomposition of width k (Theo- 
rem [63]). The algorithm can be implemented in parallel using depth polylogarithmic 
in n (Section 17. 2p . Apart from evaluating the interlace polynomial, our approach can 
also be used to compute coefficients of the interlace polynomial, for example so called 
d-truncations |Cou08l Section 5] (Section l7.3p . Our approach is not via logic but via 
the G'F(2)-rank of adjacency matrices, which is specific to the interlace polynomial. 

1.2 Obstacles 

It has been noticed that the Tutte polynomial and the interlace polynomial are sim- 
ilar in some respect |ABSn4b] : Both can be defined by a recursion using a graph 
operation, both can be defined as closed sums over edge/vertex subsets involving 
some kind of rank. These similarities suggest that evaluating the interlace polyno- 
mial using tree decompositions might work completely analogously to the respective 
approaches for the Tutte polynomial |And98t INob98] . This is not the case because 
of the following problems. 

Andrzejak's algorithm |And98] to evaluate the Tutte polynomial uses the deletion- 
contraction recursion for the Tutte polynomial (via Negami's splitting formula |Neg87| ). 
Deletion and contraction of an edge has the nice property that it is compliant with 
tree decompositions: If we are given the tree decomposition of a graph and we delete 
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Figure 1: Edge pivoting, a central graph operation for the interlace polynomial, 
increases treewidth. 

(or contract) an edge, the original tree decomposition (or, in the case of edge contrac- 
tion, a simple modification of it) is a tree decomposition of the modified graph. For 
the interlace polynomial, on the other hand, the respective graph operation is not 
compliant with tree decompositions: If we perform the pivot operation from Arratia 
et al. |ABS04a] on a graph, it is not clear how to obtain a tree decomposition of the 
modified graph. In particular, a single pivot operation can turn a tree (treewidth 1) 
into a cycle (treewidth 2), see Fig. [H 

Another problem is that in the Tutte case the recursion formula naturally gener- 
alizes from the simplest versions (chromatic polynomial) to the most general ones (it 
is the defining recursion of the Bollobas-Riordan graph invariant [BR99J : cf. also the 
recurrence relation of the polynomial of Averbouch, Godlin and Makowsky, which 
generalizes the Tutte polynomial and the matching polynomial |AGM08] ). The in- 
terlace polynomial, in contrast, needs more and more complicated recursions when 
generalizing the vertex-nullity interlace polynomial to the multivariate interlace poly- 
nomial (see Courcelle [CouOSl Proposition 12]). 

When we consider Noble's algorithm [Nob98] and concentrate on the definition of 
the Tutte / interlace polynomial by sums involving ranks, another problem emerges. 
In the Tutte case, the rank is an easy to understand graph theoretic value, namely 
the number of vertices minus the number of connected components. Noble observes 
that, if a graph is extended by a set of vertices and some edges between the old 
and the new vertices, the set of all partitions of the new vertices captures all possible 
types of "behavior" of the rank (i.e. number of connected components) when the new 
vertices and some or all of the new edges are added. - For the interlace polynomial 
on the other hand, the rank used in the definition is the rank over GF{2) of the 
adjacency matrix. Even though there exists a graph theoretic interpretation of this 
rank |Tra09] . it is substantially more involved. Furthermore, an appropriate tool to 
capture the "rank behavior" when extending a graph (such as vertex partitions in 

^But note that Traldi reduced a three-term recursion to a two-term recursion jTralO| Corol- 
lary 2.4]. 
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the case of the Tutte polynomial) seems to be missing. The main contribution of this 
work is to devise such a tool and to prove that it works well with tree decompositions. 

1.3 Outline 

We compute the interlace polynomial by dynamic programming on the tree decom- 
position of a graph. To this end, we analyze the behavior of the G'F(2)-rank of the 
adjacency matrix of a graph when the graph is extended by a fixed number of vertices 
and edges between these new vertices and the existing ones. 

Section [2] contains the definition of Courcelle's multivariate interlace polynomial, 
which we will consider in this work. We will also fix our notation for tree decom- 
positions there. In Section E] we present our approach in detail. This includes the 
motivation and definition of two central terms: extended graphs and scenarios. A 
scenario captures the behavior of the rank of an adjacency matrix when adding 
vertices. To define this precisely, we introduce symmetric Gaussian elimination in 
Section HI In Section [5l we collect properties of scenarios which enable us to use 
scenarios with tree decompositions. In Section El we describe and analyze our algo- 
rithm, which evaluates the interlace polynomial by splitting it into parts according to 
scenarios. In Section [7] we discuss how our algorithm can be parallelized and used to 
compute (some of the) coefficients of the interlace polynomial. Finally, in Section [HI 
we mention directions for further research. 

2 Preliminaries 

We consider undirected graphs without multiple edges but with self loops allowed. 
Let G = (V, E) be such a graph and A CV. By G[A] we denote the subgraph of G 
induced by A, i.e. {A, {e | e G -E, e C A}). GVA denotes the graph G with self loops 
in A toggled, i.e. the graph obtained from G by performing the following operation 
for each vertex a E A: if a has a self loop, remove it; if a does not have a self loop, 
add one. 

The adjacency matrix of G is a symmetric square matrix with entries from {0, 1}. 
As the matrices that we will consider are adjacency matrices of graphs, we use 
vertices as column/row indices. Thus, the adjacency matrix of G is a x matrix 
M = {rriuv) over {0, 1} with rriuv = 1 iSuv E E. Furthermore, we will refer to entries 
and submatrices by specifying first the rows and then the columns: the {u, f )-entry 
of M = {rriuy) is rriuv, the A x B submatrix of M is the submatrix of the entries of 
M with row index in A and column index in B. All matrix ranks will be ranks over 
the field with two elements, {0, 1} = GF(2), i.e. + is XOR and ■ is AND. Slightly 
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abusing notation we write Tk{G) for the rank of the adjacency matrix of the graph 
G. The nulhty (or co-rank) of an n x n matrix M is n(M) = n — rk(M). If G is a 
graph, we write n(G') for the nulhty of the adjacency matrix of G. 

Graph polynomials are, from a formal perspective, mappings of graphs to poly- 
nomials that respect graph isomorphism. We will consider a multivariate graph 
polynomial, the multivariate interlace polynomial. To define such a polynomial, one 
has to distinguish "ordinary" indeterminates from G-indexed indeterminates. For 
instance, x being a G-indexed indeterminate means that for each vertex a of G there 
is a different indeterminate Xa- IfACV, we write xa for HaeA^a- 

Definition 2.1 (Courcelle |Cou08j ). Let G = {V,E) be an undirected graph. The 
multivariate interlace polynomial is defined as 



where u, v are called ordinary indeterminates and x, y G-indexed indeterminates. 
2.1 Tree Decompositions 

We borrow most of our notation from Bodlaender and Koster |BK08j . A tree decom- 
position of a graph G = (V, E) is a pair ({Xj | i G /}, T = (/, F)) where T is a tree 
and each node i E I has a subset of vertices Xi (1 V associated to it, called the bag 
of i, such that the following holds: 

1. Each vertex belongs to at least one bag, that is [J^^^Xi = V. 

2. Each edge is represented by at least one bag, i.e. for all e = vw G E there is 
an i G / with v,w E Xj. 

3. For all vertices v G V, the set of nodes {i G / | f G X,} induces a connected 
subgraph of T. 

The width of a tree decomposition ({Xj}, T) is max{|Xj| | i G /} — 1. The treewidth 
of a graph G, tw(G), is the minimum width over all tree decompositions of G. 

Computing the treewidth of a graph is NP-complete. But given a graph with n 
vertices, we can compute a tree decomposition of width k (or detect that none exists) 
using Bodlaender's algorithm in time 2'^^^ |Bod96] (cf. also Downey and Fellows 




rk((GVB)[AuB]) n((GVB)[AuB]) 



(1) 



A,BQV 
AnB=0 



jDF99l Section 6.3]). 



6 



To evaluate the interlace polynomial we will use nice tree decompositions. Note 
that our definition slightly deviates from the usual on^. This has no substantial 
infiuence on the running time of the algorithms discussed in this work, but it simplifies 
the presentation. In a nice tree decomposition ({Xj},T), T is a rooted tree with 
\Xr\ = for the root r of T, and each node i of T is of one of the following types: 

• Leaf: node z is a leaf of T and =0. 

• Join: node i has exactly two children ji and j2, and Xi = Xj^ = Xj^. 

• Introduce: node i has exactly one child j, and there is a vertex a & V\Xj with 
Xi = Xj U {a}. 

• Forget: node i has exactly one child j, and there is a vertex a & V \ X^ with 
Xj = XiU {a}. 

A tree decomposition of width k with n nodes can be converted into a nice tree 
decomposition of width k with 0{n) nodes in time 0{n) ■ poly(fc) |Klo94t Lemma 
13.1.2, 13.1.3]. 

For a graph G with a nice tree decomposition ({Xj},T), we define 

Vi = ([_J{Xj I j is in the subtree of T with root i}j \ Xj and Gi = G[Vi\. 

We can think of Gi as the subgraph of G induced by all vertices that have already 
been forgotten below node i. 

3 Idea 

We will now sketch our idea how to evaluate the interlace polynomial. Our approach 
is dynamic programming similar to the work of Noble |Nob98] . Let G be a graph 
for which we want to evaluate the interlace polynomial and ({Xj},T) a nice tree 
decomposition of G. For each node i of the tree decomposition, we have defined the 
graph Gi that consists of all vertices in the bags below i that are not in Xi. We 
will compute "parts" of the interlace polynomial of Gi. These parts are essentially 
defined by the answer to the following question: How does the rank of the adjacency 
matrix of some subgraph of Gi increase when we add (some or all) vertices of Xj? 
For the leaves these parts are trivial. Our algorithm traverses the tree decomposition 

■^Usually, there is no special restriction on the bag size of the root node, and the leaf nodes 
contain exactly one vertex. 
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Figure 2: Interlace polynomial and rank behavior: Rank over GF{2) of the adjacency 
matrix increases by 2 (from 2 to 4). 



Figure 3: Interlace polynomial and rank behavior: Rank over GF{2) of the adjacency 
matrix increases by 4 (from 2 to 6). 



bottom-up. We will show how to compute the parts of an introduce, forget, or join 
node from the parts of its child node (children nodes, resp.). At the root node, there 
is only one part left. This part is the interlace polynomial of G. 

Before we go into details, let us remark that the answer to the above question 
("How does the rank of the adjacency matrix increase when adding some vertices?") 
depends on the internal structure of the graph being extended. Consider the graph 
on the left hand side in Figure [2l If we extend it by the black vertices, the rank 
increases by 2. But if we use the graph on the left hand side in Figure [Sj the same 
extension causes a rank increase by 4. 

Let us see how we handle this issue. We start with the following definition. 

Definition 3.1 (Extended graph). Let G = {V,E) be some graph, V',U C V, 
V' nU = ^. Then we define G[V',U] to denote G[V' U U] and call G[V',U] an 
extended graph, the graph obtained by extending G[V'] by U according to G. We 
call U the extension of G[V^', U]. 

Let us fix an extension U. We consider all V C V{G) such that G[V'] may be 
extended by U according to the input graph G. For any such extended graph we 
ask: "How does the rank of G[V'] increase when adding some vertices of f/?". Our 
key observation is that the answer to this question can be given without inspecting 
the actual G if we are provided with a compact description (of size independent of 
n = \V{G)\), which we call the scenario of G[V', U]. 

The scenario of (^[V^', U] (Definition 14. 4p will be constructed in the following 
way. Consider M, the adjacency matrix of G[V' U U]. Perform symmetric Gaussian 
elimination on M using only the vertices in V (for the details see Section H]). The 
resulting matrix M' is symmetric again and has the same rank as M. Furthermore, 
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Figure 4: Adjacency matrix of GfVUt/] after symmetric Gaussian elimination using 
v. Empty entries are zero. 



M' is of a form as in Figure HI The V x V' submatrix is a symmetric permutation 
matrix with some additional zero columns/rows. The nonzero entries correspond to 
edges or self loops (not of the original graph G but of some modified graph that is 
obtained from G in a well-defined way) "ruling" over their respective columns/rows: 
The edge between vi and vg rules over columns and rows vi and vg. Here, "to rule" 
means that the only Is in these columns and rows are the Is at (f i, fg) and {vs, vi). 
Similarly, the self loop at vertex rules over column and row v^. The columns and 
rows that are ruled by some edge or self loop in V are also empty (i.e. entirely zero) 
in the U xV submatrix of M'. Some columns/rows are not ruled by any edge or self 
loop in V, for instance column/row v^. This is because there is neither a self loop 
at vertex nor does it have a neighbor in V. However, may have neighbors in 
U. Thus, column of the U x V' submatrix may be any value from {0, l}*^, which 
is indicated by the question marks. Also, the contents of the U xU submatrix is not 
known to us. 

Let us choose a basis of the subspace spanned by the nonzero columns of the 
U X V submatrix and call it s^^^ . Let s^^'^ be contents of the U x U submatrix. 
By this construction, we are able to describe the rank of M' as the rank of its V x V 
submatrix plus a value that can be computed solely from s^^^ and s^^^ . 
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This will solve our problem that the rank increase depends on the internal 
structure of the graph G[V'] being extended: all we need to know is the scenario 
s = {s'^^^ ,s^^^) of (^[T^', f/]. From s, without considering (^[T^'], we can compute 
in time poly(|[/|) how the rank of the adjacency matrix of G[V'] increases when we 
add some vertices from U. This motivates the following definition. 

Definition 3.2 (Scenario). Let U be an extension, i.e. a finite set of vertices. A 
scenario of f/ is a tuple s = (s^^^ , s^^^) where s^^'^ is an ordered set of linearly 
independent vectors spanning a suhspace o/{0, 1}^ and s^^^ is a symmetric {Ux U)- 
matrix with entries from {0, 1}. A scenario for k vertices is a scenario of some vertex 
set U with \U\ = k. 

Let us come back to the evaluation of the interlace polynomial of G using a 
tree decomposition. Recall that at a node i of the tree decomposition we want to 
compute "parts" of the interlace polynomial of Essentially every scenario 

s of Xi will define such a part: The interlace polynomial itself is a sum over all 
induced subgraphs with self loops toggled for some vertices. The part of the interlace 
polynomial corresponding to scenario s will be the respective sum not over all these 
graphs but only over the ones such that s is the scenario of G[Vi, Xi]. This will lead us 
to ([2]) in SectionO To compute the parts of a join, forget and introduce node from the 
parts of its children nodes (child node, resp.), we will employ Lemma |6TT| 16. 21 and 16. 31 
These are based on the fact that scenarios are compliant with tree decompositions, 
which we will prove in Section |5] (Lemma 15. ![ Lemma 15.31 and Lemma 15. 5p . An 
example for the overall procedure of the algorithm is depicted in Figure [H 

The time bound of our algorithm stems from the following observation: The 
number of parts managed at a node i of the tree decomposition is essentially bounded 
by the number of scenarios of its bag X^. This number is independent of the size of G 
and single exponential in the bag size (and thus single exponential in the treewidth 
of G): 

Lemma 3.3. Let U be an extension, i.e. a finite set of vertices, \U\ = k. There are 
less than 2(3'^+^)*^/^ scenarios of U. 

Proof. The number of symmetric {0, l}-matrices of dimension k x k is 2^'^+^)'^/^, as 
a symmetric matrix is determined by its left lower half. Thus, there are 

2(fc+l)fc/2 

possibilities for s^^^ . 

For s^^^ , there less than 2^ possibilities: As there are 2^ — 1 non-zero elements 
of {0, 1}^, the number of linearly independent subsets of {0, 1}^ with d elements is 
bounded by ^^). Thus, the number of all linearly independent subsets of {0, l}'^ 
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Figure 5: Computation of the interlace polynomial C{G;y = 0,f = 1) of a triangle. 
In order to simplify the illustration, we ignore parameter D in which handles 
the "self loop toggling feature" of the interlace polynomial. 



is at most 



E 

0<ci<fc 



4 Symmetric Gaussian Elimination 

We want to convert adjacency matrices into matrices of a form as in Figure H] without 
touching the rank. In order to achieve this, we introduce a special way of performing 
Gaussian elimination that differs from standard Gaussian elimination in the following 
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way. First, it is symmetric, as in general every column operation is followed by a 
corresponding row operation. In this way, we maintain the correspondence between 
rows/columns of the matrix we are manipulating and vertices of a graph. Second, 
we adhere to a particular order when deciding which entry to use for the next pivot 
operation. This order is (partially) fixed by the tree decomposition. It is crucial 
for our proofs of the statements in Sect. [5] that the elimination process proceeds 
according to this order. Third, we perform symmetric Gaussian elimination using 
only vertices in a subset V of the vertices: When seeking a pivot entry in a particular 
row/column, we do not consider all entries of the row/column but only the ones that 
correspond to edges between vertices in V. 

4.1 Elimination with a Single Vertex 

Assume we are given a graph G = {V,E), its adjacency matrix M and a vertex v. 
We would like to compute the rank of M as the "effect of v on the rank" plus the 
rank of a submatrix in which we have deleted v. This might not immediately be 
possible using M itself, but we can achieve it by an appropriate modification of M. 
Arratia et al. observe that edge pivot and local complementation are such appropriate 
modifications |ABS04b] Lemma 2, Lemma 5]. For our purposes, we want to control 
the order of the operations on the adjacency matrix. Thus, we do not use edge pivot 
and local complementation directly, but define a symmetric Gaussian elimination 
step on M using v in the following way: 

• If f is an isolated vertex without a self loop, we have situation (1) of Figure El 
Vertex v has no influence on the rank of the adjacency matrix and we can 
delete the column and row corresponding to v without changing the rank of 
the adjacency matrix. The result of the elimination step is just M. 

• If V has a self loop, there is a 1 in the (f ,f)-entry of M. The elimination step 
consists of the following operations. Except for entry {v, v), we remove all Is in 
the f-column and v-row using the following pair of operations for each neighbor 
u of v: First, add the v-column to the u-column. Then, in the modified matrix, 
add the f-row to the u-row. We denote the result of the whole process by M xi f , 
which is depicted as (2) in Figure El Note that M xi f is symmetric again. The 
rank of M equals 1 plus the rank of M xi f with f-column and f-row deleted. 

Note that - up to order of the operations - this is local complementation on 
v. Writing G'" for the local complement of a graph G on vertex v |ABS04bl 
Definition 4], the adjacency matrix of is M x w |ABS04b] Proof of Lemma 
5]. 
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Figure 6: Effect of a symmetric Gaussian elimination step. Adjacency matrix with 
isolated unlooped vertex v (1), adjacency matrix after eliminating with a self loop 
at V (2), adjacency matrix after eliminating with edge vu (3). 



• If v is neither isolated nor has a self loop, there is a neighbor u of v. Assume 
that u does not have a self loop. The {u, v)- and {v, M)-entries of M equal 
1. The elimination step consists of the following operations. In the first stage, 
except for {u, v) and (f , m), we remove all Is in the f -column and f-row. This is 
accomplished by the following pair of operations for each neighbor u' of f , u' ^ 
u: First, add the -u-column to the -u'-column. Then, in the modified matrix, 
add the w-row to the w'-row. Again, performing such a pair of column/row 
operations keeps a symmetric matrix symmetric. At the end of the first stage 
the w-column and f-row consist entirely of Os, except for the entry at the 
u-position, which is 1. The second stage proceeds as follows: we add the v- 
column to every column which has a 1 in the M-row, and we also add the f-row 
to every row which has a 1 in the f-column. At the end of this stage also the 
f-column and f-row consist only of Os except at the f-position. The result of 
the second stage is a symmetric matrix again, which we denote by M xi f f . It is 
depicted as (3) in Figure [61 We do not swap columns/rows, as we must keep the 
vertices in a particular order, which is determined by the tree decomposition, 
cf. Section 14.21 The rank of M equals 2 plus the rank of M xi f f with u- and 
f -column and u- and f-row deleted. Note that this matrix is also the adjacency 
matrix of G^^\y\ {f , f }], where denotes edge pivot of G on vertices f and 
u |ABS04bi Definition 1, Lemma 2]. 

If u has a self loop we proceed analogously to obtain a matrix with a structure 
as (4) in Figure [61 Then we can eliminate the self loop at u by, say, adding 
column f to column u. (As at this point column v is zero everywhere except 
at f , only entry (u, u) of the matrix is changed by this operation and the 
symmetry is not destroyed.) Thus, we obtain a matrix exactly as (3) in Fig. [61 

We can describe the effect of a symmetric elimination step on the entries of the 
matrix (aside from the entries being set to 0) in the following way. 

Lemma 4.1. Let M = {niij) he an adjacency matrix, let a he a vertex with a self 
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Figure 7: During a symmetric Gaussian elimination step using edge ab, entry {y, x) 
is affected only by the entries at {a,x), {y,b), {y,a), {b,x) and {b,b). 



loop, and some entry of M which is not in column or row a, i.e. a ^ {x,y}. 
Then, after symmetric Gaussian elimination using a, the {y,x)-entry of M will be 

(M X a)y^ = ITLy^ + maxTriya. 

Lemma 4.2. Let M = {rriij) be an adjacency matrix, let a be a vertex without a self 
loop, ab an edge and niy^ some entry of M which is not in column or row a or b, 
i.e. {x,y} n {a,b} = 0. Then, after symmetric Gaussian elimination using ab, the 
{y,x)-entry of M will be 

(M X ab)y:r = rriy^ + niaxmyb + rUyambx + ruaxmyambb. 

We prove the statement about edge elimination, the case of self loop elimination 
is completely analogously. 

of Lemma \4.S\ Let us assume that x < y (the case x > y is analogous) . The situation 
is depicted in Figure [71 Depending on the (a,x)-entry being 1 or not, column b is 
added to column x, which adds the {y, 6)-entry to the {y, x)-entry. This gives the 
term max^riyb. After that, depending on the {y, a)-entry, row b is added to row y. 
This adds the actual value of the (6, x)-entry to the {y, x)-entry. By the previous 
column addition, the actual (6, a;)-entry is mbx + ruaxmbb- Thus, the row addition 
contributes a term mya{mbx +T^axT^bb)- The second stage has no effect on the {y.x) 
entry: Column a may be added to some other columns. But at this point of time, 
column a is entirely zero, except at the b entry. Thus, addition of the a column has 
no effect on the (|/, x) entry. The same is true for addition of the a row. □ 
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4.2 Vertex Order, Elimination with Vertex Sets, and the 
Scenario of an Extended Graph 

We want to define symmetric Gaussian elimination using a whole set V V of 
vertices. This means that we perform elimination steps using each vertex from V. 
The result of this process depends on the order in which we use the vertices for 
elimination steps. Therefore we introduce an order on the vertices of the graph, 
which will be computed before the computation of the interlace polynomial starts. 
We will use this order throughout the rest of the paper. Whenever there could be 
any ambiguity, we proceed according to this order. 

The vertex order we are using must be compliant with the tree decomposition we 
are using: Whenever a vertex is forgotten, it must be greater than all the vertices 
which have been forgotten before. Or, equivalently, the vertices in the extension Xj 
must be greater than the vertices in Vi for each node i of the tree decomposition. 
Such an order can be obtained by Algorithm [H 

Algorithm 1 Supplying a vertex order. 



1: procedure SupplyVertexOrder 

2: C 4- 1 

3: for all nodes i, in the order of bottom-up traversal, i.e. each father node is 

visited after all its children do 

4: if i is a forget node then 

5: a -h- vertex being forgotten at node i 

6: give vertex a number c in the vertex order 

7: c ^ c + 1 

8: end if 

9: end for 

10: end procedure 



Now we are ready to define elimination using a set of vertices. 

Definition 4.3. Let V V be a set of vertices of a graph G = (V, E) with adjacency 
matrix M . Symmetric Gaussian elimination on G using V is defined as the following 
process: If V = 0, we are done and M is the output of the symmetric Gaussian 
elimination process using V . Otherwise, we let v be the minimum vertex in V . If v 
has a self loop we let M' = M xi t>. Otherwise, we check whether v has a neighbor u 
in V' . If yes, we let M' = M xi vu, where u is the minimum neighbor of v. If no, we 
let M' = M. This concludes the processing of v. To complete the elimination using 
V , we continue recursively with V \ {v} in the role ofV and M' in the role of M. 
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Figure 8: Graphs corresponding to a join node i and its child nodes ji, j2- 

We also order vertex vectors (i.e. elements from {0, 1}^, U some vertex set) and 
sets of vertex vectors according to the vertex order (lexicographically). This induced 
order is used for choosing a "minimal" basis in the following definition. 

Definition 4.4 (Scenario of an extended graph). Let G\V\ U] he an extended graph 
obtained by extending G[V'] by U according to graph G = {V, E). Let the vertex order 
be such that v' < u for all v' G V and u & U. Then the scenario scen(G'[y, U]) of 
G[V', U] is defined as follows: Let M be the adjacency matrix of G[V' UU]. Perform 
symmetric Gaussian elimination on M using V to obtain M' . Let Mjjy, be the 
U X V' submatrix of M' . Consider the column space W of Mjjy,. We can choose a 
basis of W from the column vectors of M^y,. Let s^^^ be the minimal such basis. 
Let s^^^ be the contents of the U x U submatrix of M' . We define scen{G[V' , U]) 
to be (s^x^',s^><^). 

The minimal basis s^^^' in the preceding definition can by obtained by the 
following steps: Start with an empty set of columns and then as often as possible 
take the minimum column of M^jy, which is not in the span of the so far collected 
columns. 

5 Scenarios and Nice Tree Decompositions 

Consider a join node i with children ji and j2 in a nice tree decomposition of a 
graph G the interlace polynomial of which we want to evaluate. By the properties 
of tree decompositions, this implies a situation as depicted in Figure [HI Gj-^ = G[Vjj] 
and Gj2 = GlVj^] are disjoint graphs with a common extension Xj-^ = Xj^ = Xi. 
Gi = G[Vi\ = G[Vj-^ U Vj^] is the disjoint union of Gj-^ and Gj^. Assume that we have 
computed all parts (see Section [3] and ([2])) of the interlace polynomial of Gj-^ and all 
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parts of the interlace polynomial of Gj^. From this we want to compute the parts of 
the interlace polynomial of Gj. Consider one such part, say the one corresponding to 
some scenario s of X,. Somehow we have to find out for which subgraph^ G\V'] of 
Gi the scenario of the extended graph (^[T^'jXj] is s. Fortunately, these are exactly 
the subgraphs G\Vi U V2], Vi C Vj^, V2 C Vj^, with the property that the "join" of 
the scenario of G[Vi,XjJ and the scenario of G[V2,Xj2] is s. This is guaranteed by 
the following lemma. 

Lemma 5.1 (Join). Let G = {V,E) be a graph, U V , and Si,S2 two scenarios 
of U . Then there is a unique scenario S3 of U such that the following holds: If 
G[Vi] and G[V2] are disjoint subgraphs of G that may be extended by U according 
to G, scen{G[Vi,U]) = si, and seen (G [1^2, f^]) = ^2, then scen{G[Vi UV2,U]) = S3. 
Moreover, S3 can be computed from Si,S2 and G[U] within poly(|f/|) steps. 

Proof. We will apply Definition 14.41 to determine S3. We will see that S3 is uniquely 
defined by Si, S2 and G[U], and can be computed from these within the claimed time 
bound. This will prove the lemma. 

Let Gi = G[Vi] and G2 = G[V2]. Let M be the adjacency matrix of G'[ViUl^2UC/]. 
As Gi and G2 are disjoint, M has a form as depicted on the left hand side in Figure[9], 
the Vi X V2 submatrix as well as the V2 x Vi submatrix of M consists only of Os. 

By Definition 14. 4[ symmetric Gaussian elimination using Vi U V2 has to be per- 
formed on M to obtain M', which is of the form depicted on the right hand side in 
Figure [9] and from which S3 can be read off. Let us analyze a single elimination step 
occurring during the elimination process in detail, say eliminating with a self loop 
at a vertex v G Vi. One action in this step is that the 1 in the {v,v) entry will be 
used to eliminate another 1 in the v-iow by adding the v-column to the respective 
column u. Let us argue that this affects neither the Vi x V2 submatrix of M nor the 
V2 X Vi submatrix of M. As f G V^i, in the f-row the V2-entries are already 0. Thus 
we know that u ^ V2, i.e. the f-column will be added to a column from Vi U U. 
Thus, the Vi x V2 submatrix is not changed. Again as f G Vi, the V2-entries in the 
f -column are and addition of the f -column to any other column u does not change 
the V2-entries of column u. Thus, the V2 x Vi submatrix of M is not changed. 

Analogous observations can be made for the role of columns and rows reversed 
(i.e. when adding the f-row to other rows to eliminate Is in the f -column), as well as 
for elimination steps using an edge between different vertices (instead of self loops). 
We conclude that symmetric Gaussian elimination steps with Vi-vertices affect only 

^In fact induced subgraphs with self loops toggled at some vertices — but we will ignore this 
detail for the rest of the section as it is not important to understand the idea. 
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the (Vi UU) X (Vi U U) submatrix of M, but not the Vi x V2 or V2 x Vi submatrix. 
Analogously, elimination steps with V2-vertices affect only the {U U V2) x (f/ U V2) 
submatrix of M. Thus, except for the U xU submatrix, when performing symmetric 
Gaussian elimination on M using V1UV2, the same things happen as when performing 
symmetric Gaussian elimination first on G[Vi U U] using Vi and then on G[V2 U U] 
using V2. The only difference may be that depending on the vertex order elimination 
steps with Vi-vertices are interlaced with steps using V2 vertices. But we argued that 
Vi-elimination steps do not influence parts of M which are relevant for V2-elimination 
steps and vice versa, so this is not an issue. 

As elimination on M using l^i U V2 (yielding M') on the one hand does the same as 
elimination on G[Vi U U] using Vi (yielding, say, M^^^) and elimination on G[V2 U U] 
using V2 (yielding, say, M*^^)) on the other hand, the U x (Vi U V2) submatrix of 
M' is just the union of M^l^, the U xVi submatrix of Mi, and Mjjl^, the U x V2 
submatrix of M2. Recall that si^^'^^ and S2^^^^ are minimum bases of the column 
space of M^y_^, ^uv^i resp, taken from the columns of these matrices. To compute 
g^u x{ViVjV2) ^ the minimum basis of the column space of the U x [Vi U V2) submatrix 
of M' taken from the columns of this matrix, we proceed in the following way: Start 
with the empty set and as long as possible add the minimum vector of si*^^^^ □^2'^^^^ 
which is not in the span of the so far collected vectors. This can be done in time 
polynomial in \U\ using standard Gaussian elimination. 

The U xU submatrix is the only part of M which is affected by both, eliminations 
with l^i-vertices and eliminations with V2-vertices. However, the use of the U x 
U submatrix is "write-only" during the elimination process: Consider symmetric 
Gaussian elimination in general, say on some extended graph G[V' U U] using V. 
Recall that by Definition 14.31 all the elimination steps will involve only vertices from 
V in the sense that the step is either M x v ot M x vu with u, f G V. Thus, the 
contents of the U x U submatrix has no influence on what elimination steps will be 
performed. All that happens with this submatrix is that column/row vectors are 
added to it. 

Thus, the effect on the U x U submatrix of all the elimination steps during 
symmetric Gaussian elimination of G[Vi U U] using Vi can be described as adding a 
matrix, say Ai to the adjacency matrix of We can compute Ai as Ai = si^^^ — 

M{G[U]), where M{G[U]) denotes the adjacency matrix of G[U]. Analogously, we 
can compute A2 which describes the effect of symmetric Gaussian elimination of 
G'[V2U?7] using V2 on the U xU submatrix. Because of the "write-only" property, the 
effect of symmetric Gaussian elimination of M using Vi U V2 on the U xU submatrix 
of M can be described by A^ + A2. Thus we have 53^^^ = M{G[U]) + Ai + A2, 
which is the second component of S3. □ 
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Figure 10: Joining the extended graphs G[Vi, U] and G[V2, U]. 



Figure 11: Adding a vertex to an extension. 

Definition 5.2. In the situation of Lemma \5.1\ we write Sjoin(si, S2, Gff/]) for S3. 

To handle join nodes of the tree decomposition we proved Lemma f5.lt from the 
scenario of two extended graphs U] and G'[V2, U] with a common extension U we 
can compute the scenario of the joined extended graph G\Vi U V2, U] (cf. Figure [TOj) . 
To handle also introduce and forget nodes we prove two more lemmas (cf. Figure [TT| 
Figure [121). 



Lemma 5.3 (Introduce vertex). Let G = {V,E) be a graph, U ^ V , s a scenario 
of U, u & V \ U. Then there is a unique scenario s of U = U U {u} such that the 
following holds: If G[V'] may be extended by U according to G, u is not connected 
to V in G, and scen(G[V"', f/]) = s, then scen(G'[V^', [/]) = s. Moreover, s can be 
computed from s and G[U] in poly(|f/|) steps. 

Proof. As u is not connected to V',J^^^' is s^^^' with a zero component for u 
added to all the basis vectors. Also, s^^'^ is just s^^^ with a row and column added 
representing the neighbors of u in U. □ 

Definition 5.4. In the situation of Lemma \5.3\ we write s introduce {s,u,G[U]) for s. 
Figure 12: Transforming an extending vertex into a normal vertex. 
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Except for isolated vertices without self loops, every vertex has an effect on the 
rank of the adjacency matrix [ABS04bl Lemma 2, Lemma 5] (cf. Section I^?T1) . The 
following lemma states that this effect can be extracted from the scenario. 



Lemma 5.5 (Forget vertex). Let G = {V, E) be a graph, u & U ^ V , U = U \ {u}, 

V = V'U{u}, and s a scenario ofU. Then there is a unique scenario sofU and r G 
{0, 1, 2}, n G {1, 0, —1} such that the following holds: IfG[V'] is a subgraph ofG that 
may be extended by U according to G, u > v' for all v' G V , and scen(G[V^', U]) = s, 
then scen(G'[V, U]) = s and the rank (nullity) of the adjacency matrix of G[V] equals 
the rank (nullity, resp.) of the adjacency matrix of G\V'] plus r (n, resp.). Moreover, 
s and r can be computed from s and G[U] in poly(|f/|) steps, and we have n = 1 — r. 

Proof. Consider the situation after symmetric Gaussian elimination on G[V' UU] = 
G[VUU] using V (Figure [T3|) . We distinguish three cases: (1) there is a basis vector 
of the {U X V') column space with a 1 in the -u-component, (2) there is no such basis 
vector, but the (M,M)-entry of the U xU submatrix equals 1, (3) neither case (1) nor 
(2). 

Let us first consider cases (2) and (3). As all w-components of the vectors in 
s^^^' are zero, we know that symmetric Gaussian elimination on G[V U U] using 

V will consist of the following two stages: first, exactly the same operations will be 
performed as in symmetric Gaussian elimination on G[V' U U] using V' (which will 
end up in the situations depicted in Figure [13] (2), (3)), and then elimination using 
vertex u will be performed if possible. 

Thus, in case (3), s can be obtained from s in the following way: remove the 
u component of each vector of s^^^' to gain s^^^. Let a be the first column of 
s^^^. Remove the first component of a._With standard Gaussian elimination, check 
in time poly(|f/|) if a is in the span of s^^^ . If it is, let s^^^ = s^^^ ^ otherwise let 
gc/xv ^ ^f^x^ u {a}. Let s^^^ be s^^^ with first column and first row deleted. We 
have r = and n = 1. 

In case (2), we first perform an elimination step with the 1 at the (m, M)-entry: 
let s^^^ = gUxu ^ ^_ Then we continue as in case (3) but with s^^^ in the role of 
gUxu _ have r = 1 and n = 0. 

The rest of this proof deals with case (1). Let w & V he the vertex corresponding 
to the minimum vector of s^^^ with a 1 in the -u-component (cf. Figure [T^ (1)). 
Compare symmetric Gaussian elimination on ^[yut/] using V (which is performed 
to obtain s) to symmetric Gaussian elimination on G[V U U] using V (which is 
performed to obtain s). Before these two processes reach w, they are equal, but from 
w on they will differ: Using V, the edge uw will not be used for elimination and 
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Figure 13: Cases when "forgetting" an extension vertex u. Entries with question 
marks are either or 1. Empty entries are 0. 
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the process will continue with the next vertex in V' immediately. Using V, the edge 
uw will be used for elimination (which will not affect the V x V submatrix, but 
possibly change the contents of the U x (V U U) and the {V UU) x U submatrices) . 
Only after that, the process will continue with the next vertex in V. However, we 
will prove in Lemma 15.81 that we can defer the elimination using edge uw until all 
vertices of V have been proceeded and still obtain s. Thus, s can be computed in 
the following way: perform the same steps as with sjTiimetric Gaussian elimination 
on ^[K'Uf/] using V. Then, simulate the effect of a symmetric Gaussian elimination 
step using edge uw in a similar way as in cases (2) and (3). 

This simulation can be done as follows: Let w be the minimum vector of s^^^ 
with the u-component equal to 1. Let s^^^ = s^^^ \ {w} and s^^'^ = s^^^ . For 
each row i ^ u, with the Wi = 1 simulate addition of column/row u to column/row 
i doing the following: 

1. For each vector c of s^^^\ add component m of c to component i of c. 

2. Change s^^^ by first adding the u column to the i column and then, in the 
modified matrix, the u row to the i row. 

We have s^^^ = s^^^', and s^^^ is g^^^ with first column and first row removed. 
Note that after an elimination step using edge wu, the u column/row will consist 
entirely of zeros (except at {u,w) and (w,m)). Thus, the first column of s^^^ will be 
zero after the elimination with wu and we do not need to incorporate it into s^^^. 
Finally note that we have r = 2 and n = — 1 in case (1). 

□ 

Definition 5.6. In the situation of Lemma 15.51 we write Siotget{SiU,G[U]) for s, 
Arforget(s, u, G\U]) for r, and Anforgot(s, u, G[U]) for n. 

The operation defined in Definition 15.61 deletes a vertex u from a scenario in the 
sense that u is deleted from the extension but added to the graph being extended. 
We also need a notation for deleting a vertex completely from a scenario, i. e. ignoring 
some vertex of the extension. 

Definition 5.7. Let s = {s^^^\ s^^^) he a scenario of an extension U and u E U. 
Then SignoTe{s,u) is the scenario obtained from s in the following way: Delete the u- 
components from the elements of s^^^ to obtain Si. Choose the minimum (according 
to the vertex order) basis s[ for the span of si from the elements of Si using standard 
Gaussian elimination. Delete the u-column and u-row from s^^^ to obtain S2- We 
define Signore{s,u) = {s[,S2). 
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Figure 14: Symmetric Gaussian elimination using V (including steps such as elimi- 
nating with edge ah) and eliminating with edge wu can be swapped without changing 
the result. Empty entries and entries with a question mark are either or 1. 

The following lemma is used in the proof of Lemma 15.51 

Lemma 5.8. Let G = {V, E) he a graph, u E U V and G' = G[V'] a suhgraph of 
G which may he extended hy U and u > v' for all v' G V . Let w he the minimum 
vertex ofV and assume that u is the minimum neighbor of w (which implies that w 
has no neighbor in V ). Let V" = V U {u}, V = V \ {w} and M he the adjacency 
matrix of G\V' U U] (cf. Figure\T^ . Then the following two sequences of operations 
on M lead to the same result: 

1. Symmetric Gaussian elimination on M using V" , i.e. first the elimination step 
using edge wu and then the elimination steps using V . 

2. Symmetric Gaussian elimination on M using V' (i.e. the elimination steps us- 
ing V , as w has no neighbor in V ) and after that, on the result, the elimination 
step using edge wu. 

Proof. Elimination with edge wu will add the u column (row, resp.) to all columns 
(rows, resp.) which have a 1 in the w-row (column, resp.), and will then eliminate 
any remaining 1 in the u column (row, resp.). As the V-part of the w row (column., 
resp.) is entirely zero, this has no influence on the V x V submatrix of M. Thus, 
the only difference between 1. and 2. is whether the elimination step using edge wu 
is performed before or after symmetric Gaussian elimination using V. Also, it is 
enough to consider the ^/-columns and [/-rows of M. We will ignore the V x V 
submatrix of M in the following. 
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We will prove the following: every elimination step using an edge ab (a self loop at 
a, resp.) in V can be swapped with elimination using wu, i.e. the results of xia6 x wu 
and xiwu xi a6 (xa xi wu and xiwu x a, resp.) are equal. Applying this observation 
repeatedly proves the lemma. We only prove the case of an edge ab in V, the case 
of a self loop at a in ^ can be dealt with similarly. 

Let ab an edge in V. First, let us consider the column and rows of a, b, w and 
u. It is not hard to see, that, no matter whether we use first ab for elimination and 
then wu or vice versa, in the end these columns will consist entirely of zeros, except 
for {u,w), {w,u), {a,b), {b,a). Thus, it is sufficient to examine the effect of both 
elimination steps on entries (?/, x) with {x, y} H {a, 6, m, w} = 0, cf. Figure [HI 

Let M"^ = M X a6 be M after the elimination step using edge ab. Analogously 
we let M"'" = Mxiwu, as well as M"''-^'' = Mxiabxwu and M'""-"^ = Mx^wuxiab. 
We use small m to denote the entries of these matrices. For instance, m"?,'"'" denotes 
the entry in row y and column x of M"^'""". 

Case "a6 first" . By Lemma 14.21 we have 

^fx = ^yx + max ■ myb + ruya ■ rubx + rUya ■ rUax ■ rribb- 
By Lemma [4.21 again, the final value of entry {y,x) is 

ab,wu _ ab , ab _ ab , ab _ ab , ab _ ab _ ab 
'"'yx '"'yx "T '"'wx '"'yu '"'yw '"'ux '"wx '"yw '"uu^ 

where m'^^ = m^x and niy^ = my^^ as the elimination using edge ab does not 
affect column/row w (cf. Figure [T^ . Furthermore, we have 

^fu = ^yu + rriau ■ rriyb + rriya ■ nibu + man ■ niya ■ mbb, 
mfx = ^ux + max ■ mub + mua ■ mbx + max ■ mua ■ mbb, 
mtu = muu + m,au ■ m,ub + m,ua ■ m,bu + m,au ■ m,ua ■ m^bb, 

once more by Lemma [4.21 

Case "wM first". Here we have 

^«™,a6 _ ^wu _|_ ^wu _ ^wu _|_ ^wu _ ^wu _|_ ^wu _ ^wu _ ^^^"^ 

where m^" = nibb, as the entry (6, b) is not affected by edge elimination using 
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edge wu. For the remaining values we have by Lemma 14.21 
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An easy calculation yields that rriy^' = rriy^'^'^, which completes the proof. □ 

6 The Algorithm 

Algorithm [2] evaluates the interlace polynomial using a tree decomposition. The 
input for the algorithm is G = {V,E), the graph of which we want to evaluate the 
interlace polynomial, and a nice tree decomposition {{Xi}i, {I,F)) of G with 0{n) 
nodes, n = \V\. In Section [2?T] we discussed how to obtain a nice tree decomposition. 
Let k — 1 he the width of the tree decomposition, i.e. k is the maximum bag size. 

6.1 Interlace Polynomial Parts 

Algorithm [2] essentially traverses the tree decomposition bottom-up and computes 
parts S{i, D, s) of the interlace polynomial for each node i. For a node i, D O Xj, 
and a scenario s of Xj, one such part is defined in the following way: 

A,B 

where the summation extends over all A, B <Z Vi with An B = ^ and 

scen{G'[AUB,Xi]) = s, G' = GV{BUD). 

Recall that Vi is the set of vertices which have been forgotten below node i. Thus, 
S{i, D, s) is the part of the interlace polynomial of G[Vi\ corresponding to D and s. 

For every leaf i of the tree decomposition we have = and also Xi = 0. 
Thus, in Line [5] of Algorithm [2] we have D = Trivially, scen(G[0,0]) is the empty 
scenario. Thus, we have S{i, 0, ((),())) = 1 if i is a leaf. 

At the root node r the bag Xr is empty and all vertices have been forgotten, 
i.e. Vr = V . There is only one part left, S'(r, 0, ((), ()), and this is just the interlace 
polynomial of G. 
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Algorithm 2 Evaluating the interlace polynomial using a tree decomposition. 
Input: Graph G, nice tree decomposition {{Xi}i, {I,F)) of G, k such that any bag 
Xi of the tree decomposition contains at most k vertices 
1: SupplyVertexOrder > Algorithm [1] 

2: for all nodes i of the tree decomposition, in the order they appear in bottom-up 
traversal do 

for all D CXido 
if i is a leaf then 

^(z,A((),()))^i 

else if i is a join node then 

JoiN(i, D) 
else if i is an introduce node then 

lNTRODUCE(i, D) 

else if i is a forget node then 

FORGET(i, D) 
end if 
end for 
end for 

return 5(root, 0, ((), ())) > X.^ot = 
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6.2 Join Nodes 



Join nodes are handled by Algorithm [3l The correctness follows from 

Lemma 6.1. Let i be a join node with children ji and j2, D (1 Xi and s a scenario 
of Xi. Then 



where the summation extends over all scenarios Si,S2 of X^ such that 

Sjom{Sl,S2,GVD[Xi]) = s. 

Proof. Recall ([2]) for node i. Every admissible A,B give rise to Ai = A (1 Vj-^, 
A2 = AnVj^, Bi = B n Vj^, B2 = B n Vj^. G'[A U B] is the disjoint union of 
G'[Ai U Bi] and G'[A2 U B2]. (These graphs are subgraphs of the ones depicted in 
Figure El) 

We can apply Lemma [5.11 with G' in the role of G, Ai U Bi in the role of Vi and 
A2 U B2 in the role of V2. This implies that AU B takes the role of V1UV2. Using 
this it is not hard to argue that every admissible {A, B) in ([2]) corresponds to one 
pair ((v4i, Bi), {A2, B2)) of the expanded version of (I3l). □ 



Algorithm 3 Computing the parts at a join node. 
1: procedure JoiN(i, D) 
2: for all scenarios s for |Xj| vertices do 

3: > i.e., enumerate all pairs s = {s^'^^ , s^'^^^) with s^^^'^ being a list of 

linearly independent 

4: vectors from {0, 1}"^' and s^'^^' a symmetric X^ x Xj matrix with 

entries from {0, 1} 




(3) 



6 



5 



- cf. Definition 13.21 
S{i,D,s) ^ 
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8 



7: 



end for 

(ji!j2) ^ (left child of i, right child of i) 
for all scenarios si,S2 for \Xi\ vertices do 
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S ^ Sjoin{Si,S2,GVD[Xi]) 



> Definition [5]2] 
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Sii, D, s) ^ S{z, D, s) + S(ji, D, s^) ■ S{j2, D, S2 
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end for 
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end procedure 
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6.3 Introduce Nodes 

Introduce nodes are handled by Algorithm HI which is based on 

Lemma 6.2. Let i be an introduce node with child j and Xi = XjU{a}. Let D ^ Xi 
and s a scenario of Xi. Let D' = D \ {a}. Then one of the following cases applies: 

• If there is a scenario s' of Xj with Sintroduce(s', a, GV-DfXj]) = s, then we have 
S{t,D,s) = Sij,D',s'). 

• Otherwise, S{i,D,s) = 0. 

Proof. Assume there is some {A,B) such that scen{G'[AU B, Xi]) = s. Let s' = 
scen{G'[A U B, Xj]). By Lemma 1531 it follows s = Sintroducel-s', a, G'[Xi]). Conversely, 
Lemma 15.31 also guarantees that for all {A,B) with scen{G'[AU B,Xj]) = s' and 
Sintroduce(s', a, G'[Xi]) = s WB have scen(G"[y4 U B, Xi]) = s. □ 



Algorithm 4 Computing the parts at an introduce node. 
1: procedure Introduce(2, D) 
2: for all scenarios s for |Xj| vertices do 
3: S{i,D,s)^0 
4: end for 
5: j ^ child of i 

6: a ^ vertex being introduced in Xi 

7: for all scenarios s' for \Xj\ vertices do 

8: s ^ Sintroduce(s', a, (j'V-D[Xj]) > Definition 15.41 

9: Sii,D,s)^Sij,D\{a},s') 
10: end for 
11: end procedure 



6.4 Forget Nodes 

Finally, let us consider Algorithm^ which handles forget nodes. As Ariforgctis' , a, G') 
may be —1 in Lines [12] and [161 we have to assume v ^ 0. (The case f = is discussed 
in Section [6l6l ) Algorithm [5] is based on 



29 



Lemma 6.3. Let i be a forget node with child j and Xj = Xj U {a}. Let D C Xi, 
D' = D U {a} and s a scenario of X^. Then 



+ ^ 3,^yArWg,t(s',a,GVD[Xj])^Anfo,g,t(s',a,GVD[Js:j])^(^^-^ 
s'G5t 



(4) 



where 

Si = {s' I s' scenario of Xj with Signorel-s', a) = s}, 

S{ = {s' I s' scenario of Xj with Siorget{s' ,a,G'VD[Xj]) = s}, 

S{' = {s' I s' scenario of Xj with s forget W , a, GV D'[Xj]) = s}. 

Proof. We use ([2]) again. Let {A, B) be admissible. There are three cases: (1) 
a^AUB,{2)a&A and (3) a G B. In case (1), the term corresponding to {A, B) 
is contained in the first sum in (jlj). In case (2) we obtain the term corresponding 
to {A, B) from the second sum in (j4]), where we use Lemma [5.51 and multiply by Xa 
to represent the fact that a is in A. We also multiply by some power of u and v 
depending on the rank (nullity, resp.) difference with vs. without a in the extension. 
Case (3) is similar, but we also have to use D' instead of D as in this case a belongs 
to B and thus the self loop at a is toggled. □ 



6.5 Running Time 

We start with a nice tree decomposition with 0{n) nodes. Recall that k is the max- 
imum bag size of the tree decomposition. To obtain the vertex order (Algorithm [1]) 
0{n) ■ poly(fc) steps are sufficient. 

The running time of Algorithm [2] can be analyzed as follows. The i loop is 
executed 0{n) times, as there are 0{n) nodes in the tree decomposition. There are 
at most 2^ sets D <^ Xi for every node i. There are at most 

2(3fc+l)fc/2 

scenarios for 

k vertices (Lemma 13. 3p . The join case (Algorithm [3]) sums over pairs of scenarios 
and thus dominates the running time of the introduce (Algorithm H]) and forget 
(Algorithm E]) case. In the join case, we have to sum over at most (2(3'=+i)'=/2)2 

pairs 

(si,S2). Converting the scenarios (Line [10] of Algorithm |3l Line [8] of Algorithm HI 
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Algorithm 5 Computing the parts at a forget node. 

1: procedure Forget(2, D) 

2: for all scenarios s for \Xi\ vertices do 

3: S{i,D,s)^0 

4: end for 

5: j child of i 

6: a ^ vertex being forgotten in Xj 

7: for all scenarios s' for \Xj\ vertices do 

8: s Signore(s', a) > Definition 15.71 

9: S{i,D,s) i- S{i,D,s) + S{j,D,s') 

10: G' ^ G\/D[Xj] 

11: s ^ SfoTgct{s',a,G') > Definition 15.61 

12: S{i, D, S) ^ S{i, D, S) + a;„M^^fo^g-t(«'''^-^')^Anfo,get(«',a,G')5(^-^ 5') 

13: D' ^ DU {a} 

14: G' ^ GVD'[Xj\ 

15: S ^ Sforget(s',a,G") 

16: S{i, D, s) ^ S{i, D, s) + D', s') 

17: end for 

18: end procedure 



31 



and Lines [HI Eland [15] of Algorithm [5]) takes time polynomial in k, as we have shown 
in Section [5] Thus, the running time of Algorithm [2] is at most 

if we assume that arithmetic operations such as addition and multiplication (of num- 
bers) can be performed in one time step. The degree of the interlace polynomial is 
at most n in every variable (cf. Definition 12 .ip . This leads to the following result. 

Theorem 6.4. Let G = (V, E) be a graph with n vertices. Let a nice tree decomposi- 
tion of G with 0{n) nodes and width k be given, as well as numbers u, v, v ^ 0, and, 
for each a E V , Xa and Ua- Then AlgorithmlE evaluates the multivariate interlace 
polynomial C{G) at {{xa)a(^v-, {ya)a<av,u,v) using 2^^^ +o(k) . ^ arithmetic operations. 
If the bit length ofu,v, and Xa,ya,0' G V, is at most I, the operands occurring during 
the computation are of bit length 0{in). 

To evaluate the interlace polynomial of Arratia et al. |ABS04b] . which does not 
use self loop toggling in its definition, we do not need parameter D in and the 
D-loop in Algorithm [2 This simplifies the algorithm a bit. The running time is also 
reduced, but only by a factor < 2^^ and thus it is still 2'^'^ '^^^''^n. 

If we consider path decompositions (see, for example, Bodlaender jBod98] ) in- 
stead of tree decompositions, we have no join nodes. Thus, for graphs of bounded 
pathwidth, we get a result similar to Theorem 16.41 but with running time reduced to 

2l.5fc2+0(fc) . ^ 

6.6 Full-Rank Induced Subgraphs — The Case v = 0. 

If f = 0, the summation in ([1]) extends only over the A, BOV,Ar\B = ^, such 
that the adjacency matrix of GVB[A U B] has full rank. This sum can be evaluated 
using essentially the same techniques we have developed so far. Let us sketch briefly 
what changes have to be made. 

Consider the situation described on Page [H i.e. there is an extended graph 
G[V',U], and symmetric Gaussian elimination on G using V has been performed. 
The result is depicted in Figure [H Let S denote the columns of the U xV submatrix 
that are not "ruled" by any 1-entry of the V x V submatrix. (These columns are 
indicated by question marks in Figure [H) Then the following holds: The adjacency 
matrix of G[V' U U] has full rank only if S is linearly independent. If [/ = 0, the con- 
verse is also true for trivial reasons. Following this observation, we can modify our 
algorithm to count full-rank induced subgraphs only and thus evaluate the interlace 
polynomial at points with v = 0. 
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The first modification is to extend Definition 14.41 as follows: The scenario of an 
extended graph G[V',U] is said to have full rank if the column set S defined as above 
is linearly independent. 

Next, we replace ([2]) by 

^(^,D,.) = 5^x^y^«^^«^'^^)[^-^l), (5) 

A,B 

where the summation extends over all A, i? as in ([2]) with the additional restriction 
that the scenario of G'[A U B, Xi], G' = GV{B U D), must have full rank. 

Following the arguments in Section [5l it is possible to prove that full-rank scenar- 
ios can be used with tree decompositions in the same way as ordinary scenarios. For 
instance, the following version of Lemma [5.11 handles the join of full-rank scenarios: 

Lemma 6.5 (Join for full-rank). Let G = {V,E) be a graph, U 'OV, and Si,S2 be 
two scenarios of U . Then exactly one of the following cases applies: 

1. For all disjoint subgraphs G\Vi\ and G\y2\ of G such that 

(a) G\Vi\ and G\V2\ may be extended by U according to G, 

(b) G\Vi^U] has full-rank scenario Si, and 

(c) G\V2iU] has full-rank scenario S2, 

the scenario of G[Vi U V2, U] is Sjom{si, S2, G[U]) but it does not have full rank. 

2. For the same family of graphs as in the first case, the following holds: The 
scenario of G[Vi U V2, U] is Sjoin(si, S2, G[U]) and it has full rank. 

Moreover, during the po\y {\U\) -time computation of Sjoin{si, S2, G[U]) as described in 
the proof of Lemma \5.1\ it can be decided which of the two cases applies. We say 
that Sjoin(si, S2, G[U]) preserves full rank if the second case applies. 

In the algorithm, scenario-sums must be counted only if the scenario has full 
rank. For instance, join nodes can be handled by Algorithm [6l which is a slight 
modification of Algorithm [3l 

In this way. Theorem 16.41 can be established for the case f = as well. 



33 



Algorithm 6 Computing the full-rank parts at a join node. 
1: procedure JoiN_FULL_RANK(i, D) 
2: for all scenarios s for |Xj| vertices do 
3: S{i,D,s)<-0 
4: end for 

5: (ji)j2) <— (left child of i, right child of i) 

6: for all scenarios si,S2 for \Xi\ vertices do 

7: if Sjoin(si, S2, CVDlXi]) preserves full rank then 

8: S ^ Sjoin(Sl,S2,G'VD[Xi]) 

9: S{t, D, s) ^ S{i, D, s) + S{ju D, s^) ■ S{j2, D, s^) 

10: end if 

11: end for 
12: end procedure 



7 Variants of the Algorithm 
7.1 Evaluation vs. Computation 

The main motivation for our algorithm is evaluation of the multivariate interlace 
polynomial: We are given numerical values for the variables Xa,ya,u,v, an ra-vertex 
graph G and a nice tree decomposition of G. From this, we want to compute the 
numerical value C{G; {xa)aGV, {Vajaev ,u,v). Our algorithm solves this problem as 
described above. 

Another problem one might be interested in is the computation of the interlace 
polynomial: Given G, output a description of the polynomial C(G), which is a 
polynomial over the indeterminates {xa,ya I O' ^ V} U {u,v}. As the number of 
monomials of C{G) is exponential in n, there is no algorithm with running time 
polynomial in n that computes the multivariate interlace polynomial if we represent 
C{G) as a list of the coefficients of all the monomials. However, there are other ways 
of representing polynomials, for example arithmetic formulas and arithmetic circuits, 
which are considered in algebraic complexity theory |BCS97] . 

An arithmetic circuit is a directed graph with nodes of indegree or 2. Nodes 
with indegree are inputs and labeled by a constant or a variable. They compute 
the polynomial they are labeled with. Nodes with indegree two are labeled with plus 
or times and compute the sum (product, resp.) of their children. We say that a 
circuit computes a polynomial if it computes it at one of its nodes. 

If one accepts arithmetic circuits as a compact way to describe polynomials, 
then our algorithm actually computes the multivariate interlace polynomial: Use 
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Algorithm |2] as a procedure to create an arithmetic circuit for the polynomial C{G) 
in the following way. Start with a circuit with inputs Xa and ya for each a & V, a.s 
well as inputs for u, v, 0, and 1. For each operation of the algorithm of Section |6] 
using the "parts" S{i,D,S), add gates that implement this operation. In this way, 
the algorithm creates an arithmetic circuit C of size 2^'' +o{k)j^ th.aX computes C{G). 

In the following two subsections, we use this point of view for parallel evaluation 
and for computation of d-truncations of the multivariate interlace polynomial. 

7.2 Parallelization 

In this subsection we discuss a way to parallelize our algorithm. We do this using two 
operations on the tree decomposition: (1) removing all leaves and (2) contracting 
every path with more than one node. Our approach is not new but a variation of 
standard methods |Lei92l Section 2.6.1], |JaJ92t Section 3.3]. 

To describe the operations, we need some formalism. We use vectors a to collect 
the parts of the interlace polynomial which are computed. For each node i we 
define the vector ai = {S{i, D,s) \ D C Xj, s scenario of Xj), where the order of the 
components of the vector is fixed appropriately. We call cTj the "output" of node i. 
We call nodes with one child 1-nodes and nodes with two children 2-nodes. Nodes 
without children are leaves. Every 1-node has one input vector aj which is the output 
of its child, every 2-node has two input vectors which are the output vectors of its 
children. By definition, for leaves the input and the output is identical. 

With each 1-node i with child j we associate a matrix Ai. The computation of 
the 1-node i is = Aiaj. For an introduce node i with child j, by Lemma (6.21 we 
trivially can write cxj = AjCXj for some matrix A^. The entries of Ai are either or 1. 
Now let i be a forget node with child j. Consider (jl]). Note that in each of the three 
sums, the question, which S{j,D,s') {S{j,D',s'), resp.) are used, i. e. over which 
(_D, s') {{D', s'), resp.) is summed, can be answered considering only (^[Xj] and the 
involved scenarios. Thus, we can compute from this a matrix A^ with cxj = A^aj, 
too. The entries of A^ are 0, 1, Xaw'f^"' or yaU^v^~\ where / G {0, 1,2}. 

Consider a 2-node i with children ji and j2- The computation performed at i is 

(yi{D,s) = ^aj^{D,Si)aj^{D,S2), (6) 

where the sum is taken over the same elements as in (jH]). 

The parallel computation of the interlace polynomial works as follows. We start 
with the nice tree decomposition of the input graph with 0{n) nodes and an arith- 
metic circuit of constant depth which computes o", for all leaves i of the tree decom- 
position and Ai for all matrices associated with any node i of the tree decomposition. 
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Then we reduce the tree underlying the tree decomposition step by step. Every time 
we reduce the tree, we extend the arithmetic circuit such that the above invariant is 
preserved. 

We initiahze the arithmetic circuit as follows: We insert the constants and 1, 
u, V and for every vertex a of G we insert Xa and i/a- Then we produce all entries 
of all matrices associated with any node of the tree decomposition in parallel. This 
takes constant depth. 

We repeat the following operations on the tree decomposition until it consists 
only of one leaf: (1) contract all paths of 1-nodes and (2) remove all leaves. 

Path contraction works as follows. For a sequence ii,i2, ■ ■ ■ ,ii of 1-nodes we have 
(Tj^ = Ai^ ■ . . . ■ Ai-^aj, where aj is the input of node ii. Thus, we can substitute the 
sequence by one node which has A = Ai^ ■ . . . ■ Ai-^ associated with it and gets 
as input. The depth of computing the matrix product in parallel is B(log£). Thus 
a step contracting any number of disjoint 1-nodes paths of length < ^ increases the 
depth of the arithmetic circuit by 6(log£). 

Now we come to removal of leaves. By this we mean the following: Let L be the 
set of all leaves of the tree decomposition. Remove the elements of L distinguishing 
the following cases: (1) node i has two children ji and which are both leaves, (2) 
node i has two children ji and j^, one of which is a leaf (ji, say) whereas the other is 
not, and (3) node i has one child j which is a leaf. To handle case (1) we introduce 
a level with multiplications and a level with additions to perform (E]). This increases 
the depth by 2. In case (2) node i becomes a 1-node: The aj^{D,s) in ([6]) become 
coefficients of a new matrix A associated to i. As by the invariant, the arithmetic 
circuit already computes the <Jji{D, s), we do not need any new gates and depth is 
not increased. For case (3) we have to implement the matrix multiplication Aiaj to 
compute (jj. This increases the depth by a constant. Thus, removing all leaves in L 
increases the depth only by a constant. 

After performing all possible path contractions, the number of 1-nodes is at most 
two times the number of 2- nodes. Thus, at least 1/4 of the nodes are leaves. This 
implies that the following removal of leaves decreases the number of nodes of the 
tree decomposition by a factor of at least 1/4. Thus, after O(logn) steps the tree 
decomposition is reduced to a single leaf. In each step the depth increases by at most 
O(logn), which gives a O(log^n) bound on the depth of the constructed arithmetic 
circuit. 



36 



7.3 Computation of the Coefficients 



As discussed in Section ITiTl our algorithm can be used to create an arithmetic circuit 
C of size 23'='+o('=)n that computes C{G) for an n- vertex graph G with appropriate 
tree decomposition of width k. Now one can apply standard techniques to convert C 
into a procedure computing some of the coefficients of C{G). 

Let us elaborate this for an example, the computation of the d-truncation of 
the multivariate interlace polynomial. Courcelle defines the (i-truncation |Cou08t 
Section 5] of a multivariate polynomial as follows. The quasi-degree of a monomial 
is the number of vertices that index its indeterminates. As the G-indexed part of 
the monomials of the multivariate interlace polynomial are multilinear, the quasi- 
degree of a monomial of C{G) is the degree of its G-indexed part. For example, the 
quasi-degree of the monomial xaUbu^v^ is 1^41 + \B\. The d-truncation P{G)\d of a 
polynomial P{G) is the sum of its monomials of quasi-degree at most d. Let M. be 
a set of monomials. If 

/ = ^ a^m 

■m£M 

is a polynomial and Ai' C A^, we set 

f\M' = amm. 

As we want to use a result on fast multivariate polynomial multiplication which 
uses computation trees |BCS97[ Section 4.4] as model of computation, we also for- 
mulate our result in this model. In addition to the arithmetic operations (addition, 
multiplication, division), also comparisons are allowed in this model. Each of these 
operations is counted as one step. 

Theorem 7.1 ( |LS03t Theorem 1]). Consider polynomials over the indeterminates 
Xi, . . . ,Xn- Let d be a positive integer, and T) the monomials of degree at most d. Let 
f,g be two polynomials. Then, assuming the coefficients of f\V and g\'D are given, 
the coefficients of (/ ■ g)\'D can be computed using 

0{D{\ogDf\og{\ogD)) 

operations in the computation tree model, where D = 

Corollary 7.2. Let G be a graph with n vertices. Let a nice tree decomposition of 
G with width k and 0{n) nodes be given. Then the coefficients of all monomials of 
the d-truncation of C{G) can be computed using 

23fc2+0(l)^d(l+o(l))+0(l) 
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operations in the computation tree model. 

Note that the (i-truncation of C{G) has more than > n'^i^-'^osd./iogn) monomi- 
als. 



of Corollary \7.^ Let us fix a c? and a graph G with n vertices and treewidth k. 
We want to compute the coefficients of the ^-truncation of C{G). As discussed in 
Section I7.H there exists an arithmetic circuit C of size 2*^'+^('=)n computing C{G). 
We convert every operation f = g + h or f = g- h in C into a sequence of operations 
computing the coefficients of each monomial of f\d. In this way, we also get the 
coefficients of C{G)\d. To prove the corollary, it is sufficient to show that each 
operation is converted into at most n^('^~^°^^)^+'^W operations. 

We start with additions. We convert every addition gate f = g + h in C into the 
operations fm = gm + hm, m G A^, where Ai is an appropriate set of monomials. 
The monomials of C{G)\d are a subset of if denotes the set of all monomials 
over G-indexed variables x and y and ordinary variables u and v such that the quasi- 
degree is at most d and the degree in u and in v is at most n. We can select a 
monomial in Ai in the following way. First, choose d times either 1 or a variable 
from {Xa,ya I CL G V}. Then, choose the exponent of u and v from {0,1,..., n}. 
Thus, we have 

\M\ < i2n + l)'^{n + 1)2 = ^7) 

As we convert every addition from C into \Ai\ operations, the claimed bound of the 
corollary is fulfilled. 

Now let us consider multiplications, i.e. let f = g ■ h he a multiplication gate in C. 
We use fast multivariate polynomial multiplication for the G-indexed variables and 
the school method for the ordinary variables. To this end, we fix the u- and w-part 
of the monomial, i.e. we choose du and d^j, < du, d^ < n. We want to compute the 
coefficients of the monomials m of f with deg„(m) = du and deg^(m) = d^- Choose 
nonnegative integers du,g, du^h, dv^g, d^^^ such that d^^g + duA = du and d^^g + du^h = dy. 
Let 

V = {xAyB \A,BC ViG), \A\ + \B\ < d}. 

We can assume that we have already computed all coefficients of g := g\u'^^'^v'^^'^'D 
and h := h\u'^^-'^v'^^-'^V. (Here, an expression of the form u°'v^'D denotes the set 
{u°'v^m I m G T>}.) By Theorem 17. H we can compute all coefficients of the product 
g ■ h using 

0{\V\{log\V\)Hoglog\V\) = n'^V+T^J+ iog„ 
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, (if 1 I ^ 1 

operations, as \V\ < {2n + 1) < n ^ loan) , We do this for every choice of du,g-, 

du,hi dv,gi and d^^^. As these are at most [n + 1)^ many, this takes n'^{^^'^°^^^''^'' 

steps. Adding the resuhs monomial-wise needs at most |D|(n + 1)^ = n 

additions and yields the coefficients of f\u'^'^v'^^T>. We do this for all (n+ 1)^ choices 

of du and d^ to obtain the coefficients of all monomials of the (i-truncation of /. 

Thus, each multiplication in C is converted into n operations. This, 

again, is within the claimed bound of the corollary. □ 



8 Further Questions 

If we consider graphs of bounded cliquewidth instead of treewidth, so called k- 
expressions take the role of tree decompositions. Our concept of scenarios is tailor- 
made for tree decompositions and does not work with /c-expressions. Is there a linear 
algebra approach, possibly similar to the one we presented in this work, to compute 
the interlace polynomial using fc-expressions? 

The notion of rankwidth, which is related to cliquewidth |Oum05t IOS06] . is de- 
fined using the GF(2)-rank of some matrices derived from a graph. Furthermore, 
local complementation is studied in the context of the interlace polynomial as well 
as in the context of rankwidth [UumOSl Section 2]. Thus, it seems to be possible 
that rank decompositions support the computation of the interlace polynomial very 
nicely. We have not investigated this question in detail and leave it as a direction 
for further research. 
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